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Abstract 



For two vertices s and t in a graph G = (V,E), the next-to-shortest path is an 
st-path which length is minimum amongst all si-paths strictly longer than the 
shortest path length. In this paper we show that, when the graph is undirected 
and all edge lengths are nonnegative, the problem can be solved in linear time 
if the distances from s and t to all other vertices are given. 

m ; 6 

C*~) \ Keywords: Algorithms, Graphs, Shortest path, Time complexity, 

(N 

l/"") . Next-to-shortest path 

O 

1. Introduction 

Let G = (V, E, w) be an undirected graph with vertex set V, edge set E 
rS . and edge- length function w. We shall use n and m to stand for \V\ and \E\, 

respectively. For s, t G V , a simple st-path is a path from s to t without repeated 
vertex in the path. In this paper, a path always means a simple path. The length 
of a path is the total length of all edges in the path. An st-path is a shortest 
st-path if its length is minimum amongst all possible st-paths. The shortest 
path length from s to t is denoted by d(s, t) which is the length of their shortest 
path. A next-to-shortest st-path is an st-path which length is minimum amongst 
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those the path lengths strictly larger than d(s,t). And the next-to- shortest path 
problem is to find a next-to-shortest si-path for given G, s and t. In this paper, 
we present a linear time algorithm for solving the next-to-shortest path problem 
on graphs with nonnegative edge lengths, assuming the distances from s and t 
to all other vertices are given. 

History 

The next-to-shortest path problem was first studied by Lalgudi and Pa- 
mefthymiou in the directed version with no restriction to positive edge length 



10j. They showed that the problem is intractable for path and can be effi- 
ciently solved for walk (allowing repeated vertices). Algorithms for the problem 
on special graphs were also studied [2|, Il4| . For undirected graphs with positive 
edge lengths, the first polynomial algorithm was presented in 9] with time com- 
plexity 0(n 3 m) time. The time complexity has been improved several times 

12 . la [l8(. The currently best result is 0(m + nlogn) [18|, and recently the 
author further improved to linear time, assuming the distances from s and t to 
all other vertices are given. Hence, the positive length version of the next-to- 
shortest path problem can be solved with the same time complexity as the single 
source shortest paths problem. On the other hand, the problem becomes more 
complicated when edges of zero weight are allowed, and there is no polynomial 
time algorithm for this version before this work. 

Techniques 

An edge of zero-length is called as zero-edge and otherwise a positive edge. 
Let D be the union of all shortest si-paths. Let D be the digraph obtained 
from D by orientating all edges toward t. That is, for any directed edge (arc) 
in D, there is a shortest si-path in G passing through this edge with the same 
direction. Since a next-to-shortest path either contains an edge in E — E(D) 
or not, the problem is divided into two subproblems: the shortest detour path 
problem and the shortest zigzag path problem. The shortest detour path problem 
is to find a shortest si-path using at least one edge not in E(D) while the shortest 



zigzag path problem looks for a shortest si-path consisting of only edges in E{D) 
with at least one reverse arc of a positive length in D. Clearly, the shorter path 
found from the above two subproblems is a next-to-shortest path. 

In this paper, we solve the nonnegative length version also by solving the 
two subproblems individually. But there are some difficulties to be overcome. 
First, the digraph D is not so easy to construct as in the positive length version. 
Secondly, D is no more a DAG (directed acyclic graph) as in the positive length 
version, and therefore some properties in [18l | cannot be used. Instead D, we 
solve the two subproblems based on a relaxed digraph D + of D, in which all zero 
edges are regarded as bidirectional. The method to solve the shortest detour 
path subproblem is similar to the previous one for the positive length version, 
but a special care is taken into consideration for the zero-edges and the proofs 
are non-trivial and different from the previous ones. 

The shortest zigzag path subproblem is relatively more complicated. To 
solve this subproblem efficiently, the most important thing is to determine for a 
pair of vertices (a;, y) if there exists a simple si-path using a path from x to y as 
a backward subpath. The previous paper [18| showed a necessary and sufficient 
condition for the positive length version, but this condition no more holds when 
there are zero-edges. To overcome this difficulty, we use immediate dominators 
developed in the area of flow analysis. In addition, we define zero-component 
in D + , which are basically connected components of the subgraph induced by 
the zero-edges but any vertex and its dominators are divided into different 
components. By shrinking zero-components and orientating the remaining zero- 
edges, we construct an auxiliary DAG. With the help of the auxiliary DAG, 
we categorize a shortest zigzag path into four types and derive necessary and 
sufficient conditions individually. 

The main result of this paper is the following theorem, and its proof is given 
by Theorems 1151 and 1331 in Sections 3 and 4, respectively. 

Theorem 1. A next-to- shortest st -path of an undirected graph with nonnegative 
edge lengths can be found in linear time if the distances from s and t to all other 



vertices are given. 

Paper organization 

The paper is divided as follows. In Section 2, the preliminaries are presented. 
In addition to the notation used in this paper, in the preliminaries, we introduce 
dominators, a method of constructing D + , and zero-components. Also we show 
some basic properties in this section. In Sections 3 and 4, we discuss the shortest 
zigzag, and detour, path problems, respectively. And finally concluding remarks 
are given in Section 5. 

2. Preliminaries 

2.1. Notation and some properties 

Throughout this paper, we shall assume that G is the input graph and (s, t) 
is the pair of vertices for which a next-to-shortest path is asked. Furthermore, 
G is simple, connected and undirected, and all edge lengths are nonnegative 
integers. 

For a graph H, V(H) and E(H) denote its vertex and edge sets, respectively. 
For simplicity, sometimes we abuse the notation of a subgraph for its vertex set 
when there is no confusion from the context. A uv-path is a path from u to 
v. For vertices u and v on path P, let P[u,v] denote the subpath from u to v. 
We shall use "a icy-path" and a path P[u,v] alternatively. For a path P, we 
use P to denote the reverse path of P. For paths P[v\, V2} and Q[v2, V3}, P o Q 
denotes the path obtained by concatenating these two paths. Note that, even 
for an undirected path, we use P[u,v\ to specify the direction from u to v. For 
example, by "the hrst vertex x of P[u,v\ satisfying some property", we mean 
that x is the first vertex satisfying the property when we go from u to v along 
path P. Two paths are internally disjoint if they have no common vertex except 
their endpoints. For a path P, let w(P) — X) e e_E(p) w(e) denote the length of 
the path. Let d(u, v) denote the shortest path length from u to v in G, which 
is also called the distance from u to v. For convenience, let d s (v) = d(s,v) and 
d t (v) = d(v,t). 



To show the time complexities more precisely, we shall assume the distances 
from s and t to all other vertices are given. These distances can be found by 
solving the single source shortest paths (SSSP) problem. For general undirected 
and nonnegative edge length graphs (the most general setting of the problem 
discussed in this paper), the SSSP problem can be solved in 0(m + nlogn) 
time [3, |5j, and more efficient algorithms exist for special graphs or graphs 
with restrictions on edge lengths. A shortest path tree rooted at s can also be 
constructed in linear time if the distances from s to all others are given. 

2.2. D and D+ 

Let D + be the digraph obtained from D by orientating all positive edges 
toward i. That is, we treat all zero-edges as bidirectional even though only one 
direction of some of them can be used to form a shortest si-path. Our algorithm 
for finding a shortest zigzag path works on D + for the sake of efficiency. 

To construct D + , we have to construct D first. In the following, we show how 
to construct D and D + in linear time. Clearly, for v G V(D), d s (v) + dt(v) — 
d(s,t) always holds. Unfortunately, the condition that d s (v) + dt(v) = d(s,t) 
is not a necessary and sufficient condition to determine the set of vertices in 
V(D) when there are zero-edges. The reason is described as follows. Let D' 
be the subgraph of G with V(D') = {v\d s (v) + d t (v) = d(s,t)} and E(D') = 
{(u, v)\d s (v) = d s (u) + w{u, v)} for u, v G V(D'). A vertex is a non-st-cut if it 
is a cut vertex and its removal does not separate s and t. For a non-si-cut x, a 
connected component K of D' — x is called a knob if s,i ^ V(K). Since a; is a 
cut vertex, any si-path passing through a vertex in K repeats at x and cannot 
be simple. Furthermore, for any vertex v in K, since d s (v) + dt(v) — d(s,t), it 
must be connected to a; by a path of zero-length. 

Lemma 2. A vertex v is in V(D) iff v G V(D') is not in any knob. 

Proof. By definition, v G V(D) implies v G V(D'). Furthermore, v cannot be 
in any knob since there is no simple si-path in D' passing through any vertex 
in a knob. 



Now, we prove the other direction. For any vertex v G V(D') ~ {s,t}, con- 
sider the digraph D" obtained from D' by reversing the direction of all positive 
edges (x, y) with d s (y) > d s (v). Also we add a new vertex sq as well as two arcs 
(so, s) and (so, t). Then there exists a shortest st-path passing v in D' iff there 
are two disjoint paths from so to v in D" ', or equivalently there is no non-si-cut. 
Obviously any vertex u with d a [u) 7^ d a (u) cannot be an sof-cut in D" , and 
there exists such a cut node iff v is in a knob. D 



Lemma 3. D + can be constructed in linear time if d s (v) and dt(v) are given 
for all v. 

Proof. First we construct D' in linear time. By using depth-first search start- 
ing from s, all cut vertices can be found in linear time. According to the order 
of found cut vertices, all knobs can be detected by checking the components 
after removing the cut vertices. D 



2.3. Dominators in D + 

We shall use the term "immediate dominators" defined in [lj. A vertex 
v G V(D + ) is an s-dominator of another vertex u iff all paths from s to 11 
contain v. An s-dominator v of u is an s -immediate- dominator of u, denoted 
by I s (u), if it is the one closest to u, i.e., any other s-dominator of u is an 
s-dominator of I s (u). In D + , any vertex has a unique s-immediate-dominator. 
The i-dominator is defined symmetrically, i.e., v is a t-dominator of u iff any 
ui-path contains v, and It{u) stands for the i-dominator closest to u. Note that 
s is an s-dominator and t is a t-dominator of any other vertex in D + . 

Finding immediate dominators is one of the most fundamental problems in 
the area of global flow analysis and program optimization. The first algorithm 
for the problem was proposed in 1969 b y L owry and Medlock [13J], and then had 
been improved several times [y, I 111 . Ila . Il6| . A linear time algorithm for finding 



the immediate dominator for each vertex was given in 
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2.4- Zero- components 

Definition 1. A path P is a 0-path if all edges in P are zero-edges. A 0-path 
P[u, v] is a 0*-path if P[u, v] does not contain any vertex in {I s (u), I 3 (v), It(u), It(v)}. 
A zero- component is the subgraph of D + induced by a maximal set of vertices 
in which every two vertices are connected by a 0*-path. The zero-component 
which v belongs to is denoted by Z(v). 

A zero-component may contain only one vertex but no edge. All the zero- 
components partition V(D) into equivalence classes, i.e., v\ £ Z{vi) iff v-i e 
Z(v\). We shall show how to find all zero-components of D + in linear time. 

Lemma 4. If v' e Z(v), then I s (v') = I s (v) and /*(«') = It{v). 

Proof. If I s (v') is not an s-dominator of v, there is an su-path Q\ avoiding 
I s (v'). Since v and v' are in the same zero-component, there is a 0-path Q2[v, v'] 
in Z(v) avoiding I s (v'). Thus, Q\ o Q 2 , possibly taking a short-cut if the path 
is non-simple, is a path from s to v' avoiding I s (v'), a contradiction. Therefore 
I s {v') is an s-dominator of v. Similarly we can show that I s (v) is also an s- 
dominator of v'. Consequently I s (v') and I s (v) dominate each other, and thus 
they are the same vertex. The result It{v') — It{v) can be shown similarly. □ 

An s-dominator tree [l| of D + is a tree T with root s and vertex set V(D + ). 
A vertex u is a child of v in T iff v = I s (u) . 

Lemma 5. The subgraph of D + induced by the edge set E — E{Tf) — E(Tf) 
is the union of all zero- components, where Eq is the zero-edges set, and E{Tf) 
and E(Tf) are the edge sets of s- and t-dominator trees of D + , respectively. 

Proof. Since no positive edge is in any zero-component, we only need to con- 
sider the zero-edges Eq, For any vertex v, if (u,I s (v)) is the last edge of a 



path from v to I s (v), then u is a child of I s (v) in the s-dominator tree. Af- 
ter removing E(T^) and E(Tf), there is no path from any vertex to its s- or 
i-dominator. Therefore, by definition, the induced subgraph is the union of all 
zero-components. □ 

Since a dominator tree can be constructed in linear time |l| , the next corol- 
lary follows directly from the above lemma. 

Corollary 6. All zero- components of D + can be found in linear time. 

2.5. Outward and backward subpaths 

A positive edge (u,v) g E is a reverse positive edge if (v,u) € E(D + ). It 
implies that (u, v) ^ E(D + ) since any positive edge in D + is unidirectional. 

Definition 2. A backward subpath of a path in G is a path consisting of at least 
one reverse positive edge and possibly some zero-edges. A semi-path in D + with 
at least one backward subpath is called a zigzag path. Two backward subpaths 
in a zigzag path are consecutive if there are separated by a sequence of non- 
reverse positive edges and zero edges; otherwise, they form a longer backward 
subpath indeed. 4 

By definition, a zigzag path is a semi-path in D + . For simplicity, we shall 
use "path" instead of "semi-path" in the following. 

Definition 3. An outward subpath of an si-path in G is a path consisting of 
edges in E — E{D). The both endpoints of an outward subpath are in V(D) and 
all its internal vertices, if any, are not in V(D). An si-path is called a detour 
path if it contains at least one outward subpath. 



4 Another way to define a backward subpath is a maximal subpath consisting of at least one 
reverse positive edge and possibly some zero-edges. The difference is that, by our definition, 
there may be some zero-edges preceding or succeeding a backward subpath. Our definition is 
for the sake of simplifying some proofs. 



The shortest detour path problem is to find a shortest detour si-path while 
the shortest zigzag path problem looks for a shortest zigzag st-path consisting 
of only edges in E(D). Since a next-to-shortest path either contains an edge 
in E — E{D) or not, the shorter path found from the above two subproblems 
is a next-to-shortest path. Since s and t are fixed throughout this paper, we 
shall simply use "zigzag path" and "detour path" instead of "zigzag st-path" 
and "detour st-path" , respectively. 

When the edge lengths are all positive, the following result was shown in 
and it is also the basis of the algorithms in this paper. In remaining paragraphs 
of this subsection, we show Theorem [7] by Lemmas [5] and [TU1 

Theorem 7. A shortest zigzag path contains exactly one backward subpath. A 
shortest detour path contains exactly one outward subpath and no backward sub- 
path. 

Lemma 8. A shortest zigzag path contains exactly one backward subpath. 

Proof. Suppose by contradiction that P is a shortest zigzag path in D + with 
more than one backward subpath. Let P[xi, yt], for 1 ^ i ^ fc, be the consecutive 
backward subpaths in P and Q = P[xx, yk] where d s (Xi) > d s {yi) and k Jj 2. We 
may assume that the first and the last edges of Q are positive edges (otherwise 
move X\ forth or y^ back accordingly). Let x' be the first vertex on P such 
that w(P[x',Xi\) = and y' the last vertex such that w{P[yk,y'\) — (see 
Fig. [TJ(a)). We divide into three cases, and in either case we show that there 
exists a shorter zigzag path P'. 

• There is a path Pi from s to an internal vertex v of Q such that P\ is 
disjoint to P[yk, y']- Then P' — Pi o P[v, t] is a zigzag path. Since Pi is a 
short-cut of P[s, v], P' is shorter than P (see Fig. 03(b)). 

• There is a path P 2 from an internal vertex v of Q to t such that P 2 is 
disjoint to P[x',xi\. Similarly, P' — P[s,v] o P 2 is a zigzag path shorter 
than P. 





Figure 1: Illustrations for Lemma [8] (a) A path P with more than one backward subpath. 
The bold line is Q; (b) Case 1; (c) Case 3. 



• Otherwise, since the first case does not hold, there exists a path P\ from s 
to a vertex v\ on P[yk,y'}, which is internally disjoint to Q. Furthermore, 
d s (y') $J d s (yi) < d s (x'). Similarly, there exists a path Pi from a vertex vi 
on P[x' , xi] to t, which is internally disjoint to Q. And dt(x') ^ dt(Xk) < 
dt(y'). Then the path P' = P\ o P[v\,vi] ° Pi is a zigzag path. Clearly 
w(P') = w(P) - 4( x ') - dt(y') + d s {y') + dt(x') < w(P) (see Fig. ffl(c)). 

D 



Lemma 9. For any two vertices x and y in V(D + ) — {s,t}, there exist an 
sx-path and a yt-path; or an sy-path and an xt-path; which are disjoint. 

Proof. The result is trivial for the case x £ Z(y). We only need to show the 
case x £ Z(y). To show the lemma for this case, we construct an auxiliary 
directed graph from D + by adding a new vertex v and two bidirectional edges 
(v,x) and (v,y). Since there is no non-si-cut, similar to Menger's theorem, 
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there is an si-path passing through v in the auxiliary graph, and the desired 
two paths exist. □ 



Lemma 10. A shortest detour path contains exactly one outward subpath and 
no backward subpath. 

Proof. Let P be a shortest detour path, in which P[x, y] is an outward sub- 
path. We shall show that if P had another outward subpath or backward sub- 
path in addition to P[x,y), we could construct a detour path P' shorter than 
P. 

By Lemma [5J there exist an sx-path and a j/i-path; or an sy-path and an 
a;i-path in D + which are disjoint. In either case that the two paths exist, we can 
concatenate the two paths with P[x, y] (or its reverse) to form a simple si-path. 
It is clear that the shorter detour path in the two cases is a shortest detour path 
P' containing P[x, y]. D 



3. Shortest zigzag path 

3.1. Basic properties 

By Theorem |51 a shortest zigzag path has the form P* = Pi [s, x] o P 2 [x, y] o 
P^[y,t], in which Pi are paths in D + . Since P* is required to be simple, the 
three subpaths must be simple and disjoint except at the two joint vertices. 
Therefore our goal is to find i,y£ V(D) minimizing 

d(s, x) + d(x, y) + d(y, t) = d(s, t) + 2d(y, x) (1) 

subject to that there exists a simple path Pi[s,x] o P 2 [x,y] o P 3 [y,t] in D + . 
Since d(s,t) is fixed for a given graph G, the objective is to find the minimum 
of d(y,x). If x and y satisfy the constraint, we say "the pair (x,y) is valid^ 
and "y is valid for x" . A valid pair (x, y) with minimum d{y, x) is an optimal 
backward pair, or simply optimal pair, and the corresponding backward subpath 
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is an optimal backward subpath. The shortest zigzag path problem is equivalent 
to finding an optimal pair. 

The auxiliary simple digraph Z is obtained from D + by shrinking every 
zero-component and orientating all the remaining zero-edges toward t. By the 
definition of zero-component, if (u, v) is a zero-edge in Z, u = I s (v) or v = It{u). 
Therefore the orientation can be easily done. For Z(v) in D + . let z v denote its 
corresponding vertex in Z. For simplicity, since s and t themselves must be 
zero-components, the corresponding vertices in Z are also denoted by s and t, 
respectively. For a vertex z v , I s {z v ) and It(z v ) are again the immediate s- and 
i-dominators (but in Z). A simple path in D + corresponds to a simple path 
in Z since, without backward subpath, a path cannot enter a zero-component 
twice. 

Definition 4. We define a binary relation on pairs of vertices in V(D + ): u < v 
or equivalently v y u iff z u ^ z v and there exists a path from z u to z v in Z. Let 
C s (u) = {v\I s (u) ~< u A u -< u} and C t (u) = {w|w -< I t (u) A u -< v}. 

Definition 5. The predicate Pi(x,y) is true iff y G C s (x) and x £ C t (y). 

Lemma 11. If P\{x,y) is true, d s (y) < d s (x). 

Proof. By definition, y G C s (x), and therefore d s (y) ^ d s (x) and y ^ Z(x). If 
d s (y) = d s (x), they are connected by a 0-path but not a 0*-path, i.e., a path con- 
taining a vertex in {/ s (a;),/ s (j/),/t(a;),7t(2/)}. Since y -< x, a yx-path contains 
neither 7 8 (j/) nor Jf(x). Since y G C s (x) and x & C t {y), I s {x) < y -< x -< I t (y), 
which implies that any yx-path in Z? + contains neither / s (a;) nor h{y), a con- 
tradiction. D 

The notation defined on D + will also be used for Z. We do not distinguish 
between them since there will be no confusion from the context. The next two 
lemmas appeared in [18l | for positive length version, and it is easy to see it 
also holds for nonnegative length version. The next lemma show a necessary 
condition for the validity of a pair. 
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Lemma 12. If (x,y) is valid, then f5\(x,y) is true. 

Proof. By definition, y -< x. If I s (x) -/< y, by the definition of immediate 
dominator, any sx-path and yx-path contain I s (x) simultaneously and cannot 
be disjoint. Therefore we have I s (x) -< y, and then y £ C s (x) by definition. 
The relation x £ Ct (y) can be shown similarly. □ 

Lemma 13. If y £ C s (x), there are two paths from s and y, respectively, to x, 
which are disjoint except at x. 

Proof. Let p = I s (x) and R be any sp-path. By the definition of immediate 
dominator, removing any vertex in C s {x) cannot separate p and x and therefore 
there are two internally disjoint px-pedhs, say P\ and P 2 . If y is on one of them, 
say P2, we have done since RoP\ and P 2 [y, x] are the desired paths. Otherwise, 
let P3 be any yx-p&th and v be the first vertex on P3 and also in V(P\) U V(P2)- 
W.l.o.g. let v £ V(P\). Then, the path P%\y, v] o Pi[v,x] is a yx-p&th disjoint 
toRoP 2 . D 

Lemma 14. If (3\(x,y) is true and there exists a path P from y to t avoiding 
Z(x), then there exists a vertex valid for x. Furthermore if y* satisfies the above 
condition with minimum d(y*,x), then there exists a vertex v such that (x,v) is 
valid and d(v, x) = d(y* , x). The same result also holds for the case that /3i(x, y) 
is true and there exists a path from s to x avoiding Z(y). 

Proof. We show the first result and the second one can be shown similarly. 
Let v be the last vertex of P in C s (x). Since fli(x,v) is also true, we have that 
d s (y) ^ d s (v) < d s (x) by Lemma 1111 By Lemma 1 131 there are a path Pi[s,a;] 
and a path P 2 [v, x] which are internally disjoint. Then, the path P\o P 2 o P[v,t\ 
is a zigzag path and therefore (x, v) is a valid pair. 

Since v is also a vertex satisfying the condition, we have d(v,x) = d(y*,x), 
otherwise v contradicts the minimality of y* . □ 
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Types of optimal backward pairs 

By the definition of zero-component, there exists an sx-path avoiding Z(y) 
iff z y is not an s-dominator of z x . Similarly, there exists a yi-path avoiding 
Z(x) iff z x is not a t-dominator of z y . Therefore, all the valid pairs (x,y) can 
be categorized into the following four types, and the best of the four types, if 
any, is an optimal pair. 



• 



• 



• 



• 



Type I: z y is not an s-dominator of z x and z x is not a i-dominator of z y . 



Type II: z y is an s-dominator of z x and z x is a i-dominator of z. 



v 



Type III: z y is an s-dominator of z x and z x is not a i-dominator of z, 



v 



Type IV: z y is not an s-dominator of z x and z x is a t-dominator of z y . 



In the following subsections, we shall derive linear time algorithms for each 
of the types. The next theorem concludes the result of this section, and its proof 
is given by Lemmas [TBI [25l and [28l in the following subsections. 

Theorem 15. Suppose that d s (v) and dt(v) are given for all vertices v. A 
shortest zigzag path can be found in linear time. 

3.2. Type I 

By definition, Z is a DAG. If (z u , z v ) is a zero-edge in Z, then z u = I s (z v ) 
or z v = I t (z u ). By this prop erty , all the properties and the algorithm derived 



for a shortest zigzag path in 



181 also hold for Z. 



Lemma 16. Suppose that d s (v) and dt(v) are given for all vertices v. A short- 
est zigzag path of type I can be found in linear time. 

Proof. For (x,y) such that /?i(x,y) is true, by definition, the pair {x,y) is 
valid of type I iff {z x , z y ) is valid in Z. Therefore, a shortest zigzag path of type 
I in D + can be found by solving the shortest zigzag path problem in Z. By the 
result of [18], it can be done in linear time. 4 □ 
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(a) 



(b) 



Figure 2: Illustrations for Lemma ll7l The dotted-line depicts a valid zigzag path, (a) V(Q) n 
V(P) = and (x',q) is valid; (b) Q intersects Pi and (v,y) is valid. 



3.3. Type II 

For a shortest zigzag path of type II, the corresponding path in Z repeats 
at both z x and z y . The next lemma is not only for type II. 

Lemma 17. For any y G C s {x), if y is valid for some x' >- x, then there exists 
some v £ Z{x) such that (v,y) or (x',v) is valid. 

Proof. Since y is valid for x' , there exists a path P = Pi[s, x'jo P 2 [x' ,y]oP 3 [y,i\. 
By Lemma fl2l f3i(x\y) is true and x' -< h(y)- Since y E C s {x) and i-!a/X 
It(y), Pi{x,y) is also true. 

If P<i or P3 does not pass any vertex in Z(x), by Lemma [14l (x, y) is valid and 
the proof is complete. Otherwise both the two subpaths pass vertices in Z(x), 
and therefore, in Z{x) we can find q and q' on Pi and P3, respectively, as well 
as a 0-path Q[g, </] which is internally disjoint to Pi and P3. Since y G C s (x') 
and y -< a; -< a/, a; G C s (x'). Since P3 is a path passing Z(x) and avoiding a;', 
a;' G Ct(x). Therefore d s (x') > d s (x) by Lemma [TT1 



*The problem is named the optimal backward problem in [1 
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If Q is disjoint to Pi, the path P = Pi[s, x']oP 2 [x' ,q] oQoP 3 [q',t) is a simple 
path with a backward subpath from x' to q. That is, q € Z(x) and (x',q) is 
valid (Fig. 0(a)). Otherwise Q intersects P\. Let v be the intersection vertex 
closest to q on Q. Then, the path P — Pi[s,v] o Q[v,q] o P 2 [q,y] o P 3 [y,t] is a 
simple path with a backward subpath from v to y. That is, v € Z(x) and (v, y) 
is valid (Fig. (3(b)). □ 



Lemma 18. For any two vertices u and v such that I s {u) = I s (v) = p, there 
exist two internally disjoint paths from p to u and v, respectively. 

Proof. By definition, there exists no cut vertex whose removal separates I s (v) 
from u or v. By Menger's theorem, such two disjoint paths exist. □ 

By definition, if (x, y) is valid for type II, z x is a t-dominator of z y and z y is 
an s-dominator of z x . We show a stronger condition in the next lemma. 

Lemma 19. If(x,y) is an optimal pair of type II, z x = It(z y ) and z y — I s (z x ). 

Proof. Suppose that P = P\ o P 2 o P 3 is a shortest zigzag path of type II, 
in which P 2 is the backward subpath from x to y. If z y ^ I s (z x ) — z v , both 
Pi and P 2 contain a vertex in Z(v). We shall show that y £ C s (v), and then 
by Lemma ITTl (x,y) is not optimal. The result z x — It(z y ) can be handled 
similarly. 

Let P{[yi,vi) and P 2 [3/2,^2] be subpaths of P\ and P2, respectively, such 
that {2/1,2/2} C Z(y), {i>i,i>2} C Z(v), and no internal vertex of them is in 
Z(y) U Z(v). We can find a path Q[fi,W2] in Z(v) and two disjoint paths from 
I s (y) to j/i and 2/2, respectively (Lemma ITS)) . Then there are two disjoint paths 
from vi to I s (y), and therefore y 6 C s (v). □ 



If z x = It(z y ) and z y — I s {z x ) as well as /3i(x,y) is true, we say that Z(x) 
and Z{y) are a candidate component pair. By Lemma I19[ to find an optimal 
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Figure 3: (a) Four cases of the three paths when /32(x,y) is true; (b) The case that Ri and 
R2 exist (Lemma [22}; (c) The case that R3 and R4 exist (Xemma l22| l, 



pair of type II, we only need to determine if there exists a valid pair for any 
candidate component pair. For a candidate component pair Z(x) and Z(y), let 
H xy be the digraph with vertex set U — {v\y < v hv < x}U Z(y) U Z(x). The 
edge set is E(D+[U]) - E{Z(x)) - E(Z(y)), in which D+[U] is the subgraph of 
D + induced by U . 

Definition 6. The predicate j3 2 (x,y) is true iff there are three internally dis- 
joint paths Pi[yi, Xi], P 2 [y 2 ,x 2 ] &ndP 3 [y 3 ,x 3 ] iaH xy satisfying: (1) {x 1 ,x 2 ,x 3 } C 
Z(x) and {2/1,2/2,2/3} C Z(y); and (2) x t ^ Xj or y { ^ y 3 for 1 ^ i,j ^ 3; and 
(3) 2/1 ^ 2/2 and x 2 ^ £3- 

Fig. [21(a) illustrates the four possible cases of the three paths when f3 2 (x, y) 
is true. 

Lemma 20. For any three vertices y\, y 2 and 2/3 in Z(y), there are two disjoint 
paths P[s, 2/1] andQ[y 3 ,y 2 ], or P[s 1 y 2 ] andQ[y 3 ,yx}. 

Proof. Let p — I s (y) and R be any sp-path. Note that R avoids Z(y). By 
Lemma 1181 there are two disjoint paths P\[p,y\\ and P 2 [p, 2/2]- If 2/3 i s on one 
of the paths, we have done. Otherwise, let P 3 be any path in Z(yi) from y 2 to 
2/3, and q be the last vertex of P 3 appeared on P\ or P 2 . If q is on P\, then 
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Q = P\ [yi, g] o Ps[q, ys] and P = R o P 2 arc the desired two paths. The case 
that q is on P 2 can be shown similarly. □ 



Corollary 21. For any three vertices x\, x 2 and x% in Z{x), there are two 
disjoint paths P[x 2 ,t] and Q[x±, X3], or P[x 3 ,t] and Q [x± , X2] ■ 

Lemma 22. Suppose thatZix) andZ(y) are a candidate component pair. There 
exist x' £ Z{x) and y' £ Z{y) such that (x',y') is a valid backward pair of type 
II iff f3 2 (x, y) is true. 

Proof. It is clear that if (a/, y') is valid for type II, the three paths exist and 
fi 2 (x, y) is true. It remains to prove that such a valid (x' , y') exists if f3 2 (x, y) is 
true. Since (3 2 (x,y) is true, there are three internally disjoint paths Pi [2/1,21], 
Pi [y2,x 2 ] and P 3 [y 3 , x 3 ] in H xy , in which yi ^ y 2 and x 2 ^ x 3 . 

By Lemma [20l and w.l.o.g., there exist disjoint paths Qi[s,y{\ and Q2[j/3,2/2]. 
In the case that y 2 — j/3, Q2 contains only one vertex but no edge. If x\ ^ x 2 , by 
Corollary 12 11 there exist two disjoint paths Ri[x 2 ,t] and R 2 [xi^x 3 ]; or R 3 [x3,t] 
and -Pm^i, x 2 }. If Pi and R 2 exist, the path Q 1 o P 1 o R 2 o P 3 o Q 2 o P 2 o R x is 
a desired path, as shown in Fig. (3](b). That is, x' = x\ and y' — j/3. Otherwise 
R3 and R4 exist, and the path Qi o Pi o P 4 o P 2 o Q 2 o P 3 o P 3 is the desired 
path, as shown in Fig. [31(c), namely, x' = x\ and y' — y 2 . 

It remains to consider x\ = x 2 . By the definition of immediate dominator, 
there is a path R from X3 to t avoiding x\. The path Qi o Pi o P 2 o Q 2 o P 3 o R 
is a desired path (similar to Fig. [21(c)), in which x' ~ X\ and y' — y 2 . D 



From H xy , we construct a vertex-capacitated digraph ip+ as follows. First, 
if any connected component contains exactly one vertex u in Z(x) and one 
vertex v in Z{y), we replace the component by an edge (v,u). Then, we add a 
new source ya and a new sink xq. For each v G Z{y) there is an edge (j/Oj v); and 
there is an edge (u, xq) for each v £ Z(x). The capacity of vertex v is denoted 
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Figure 4: Illustrations for Lemma 1231 (a) Q and P3 are disjoint; (b) Q intersects Pg. The 
dotted-line depicts a valid zigzag path. 

by c(v). The capacities are assigned as follows: c(xo) = c(j/o) = °°; c i v ) = 2 
for any v £ Z[x) U Z(y); and c(u) = 1 for any other vertex. In the following, 
"the max-flow in Hj" means the maximum vertex-capacitated flow from j/o to 
xq in H^. Note that a vertex-capacitated digraph can be easily transformed 
to an edge-capacitated digraph, and the maximum flow of a vertex-capacitated 
digraph can be computed by traditional maximum-flow algorithms. 

Lemma 23. Suppose that Z(x) and Z{y) are a candidate component pair and 
(x',v) is not a valid type-I pair for any x' G Z(x) and d s (v) Js d s (y). Then, 
faix, y) is true iff the max-flow in H^ y is at least three. 

Proof. If ^(x, y) is true, it is easy to see that the max-flow in 7J+ is at least 
three. We need to show the other direction. If there is a flow of value three, 
there are three internally disjoint paths Pi[yi,Xi], 1 ^ i ^ 3. According to 
the assigned capacities, there are at least two distinct vertices in {2/1,2/2,2/3}, 
and so are in {x\,X2,X3}. The only question is that two of the three paths 
may have the same endpoints. That is, w.l.o.g., x\ = X2 and 2/1 =2/2- By the 
construction of Hj y , the connected component containing the two paths must 
contain another vertex in Z(x) U Z(y). W.l.o.g. let y' 6 Z(y) be such a vertex. 
Then, let Q be a path from y' to x± and q the first vertex of Q intersecting P± 
or P2. W.l.o.g. let q be on Pi (see Fig. SJ. 
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• If Q and P3 are disjoint, the three paths Pi, Q[y',q] o P 2 [q,Xi] and P3 
satisfy the requirement and f5<2{x,y) is true (Fig. 01(a)). 

• Otherwise Q and P3 share a common vertex, possibly y' — j/3. Let q' be 
the last vertex of Q on Pj, . 

— If q' = 2/3, the three paths Pi, Q[y3,q] o P 2 [q,xi] and P3 satisfy the 
requirement and ^ 2 (x,y) is true. 

— Otherwise d s (q') ^ d s (y). There exists a path Pi oP 2 [ii, q] °Q[q, q'] ° 
P3 [q' , X3] is a path from ?/i to X3 with a backward subpath of length 
d(q',Xi) ^ d(y,x), and this path can be extended to a zigzag path 
of type I, a contradiction to the assumption (Fig. S(b)). Note that 
(3i(x,q') is true and therefore d s {q') < d s (x) by Lemma [HI 

□ 



Corollary 24. Under the assumption of Lemma \23l ^2(^,2/) can be determined 
in 0(m xy + n xy ) time, in which m xy and n xy are the numbers of edges and 
vertices in Hj, respectively. 

Proof. By Lemma [231 P2,{x,y) can be determined by checking whether the 
max-flow in H xv is larger than or equal to three. Since all the capacities are 
integral, this max-flow question can be determined with at most three iterations 
of the augmentation step of the Ford-Fulkerson maximum flow algorithm [J, |3j 
or equivalently at most three breadth-first search on the residue graphs. There- 
fore the time complexity is linear. □ 



Lemma 25. Suppose that I is the length of an optimal backward subpath of type 
I. In linear time, we can find an optimal backward subpath of type II with length 
less than I or determine there is no such subpath. 
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Proof. By Lemmas [5U and [231 we can determine if there exists an optimal 
backward subpath of type II with length less than I. Note that the proofs of 
Lemmas [2D and [531 are constructive and they implies a linear time algorithm for 
constructing such a zigzag path if it exists. By Corollary |24l the time complexity 
is linear to ^2(m xy + n xy ), in which the summation is taken over all candidate 
component pairs. By Lemma \W\ and the uniqueness of immediate dominator, 
any zero-component is involved in the max-flow computations at most twice. 
Therefore the total time complexity is 0(m + n). □ 



3.4. Types III and IV 

Types III and IV are similar to Type II, but simpler. Furthermore, the two 
types are symmetric and we shall only explain Type III briefly. A pair (x, y) is 
valid for type III if z y is an s-dominator of z x and z x is not a i-dominator of z y . 

Lemma 26. If (x,y) is an optimal pair of type III, z y = I s (z x ). 

Proof. By using a similar argument as in Lemma \T§[ this lemma follows. □ 

In the next lemma, H xy has the same definition as in type II. 

Lemma 27. Suppose that z y = I s (z x ) and z x is not a t-dominator of z y . There 
exist x' G Z{x) and y' £ Z{y) such that (x',y') is a valid backward pair of type 
III iff there are two disjoint paths Pi[yi,xi] and P2 [2/2, X2] i n H xy such that 

Proof. It is clear that if (x',y') is valid for type III, the two paths exist. 
Conversely, if Pi and Pi exist, by Lemma Q~51 there are two disjoint paths 
Qi[I s (y),yi] and Q2[Is(y),y2], respectively. Since z x is not a t-dominator of z y , 
we can find a path R from y\ to t and avoiding Z(x). Let v be the last vertex 
of R intersecting Q\ or Q2. Then v is valid for x\ , namely, x' = x\ and y' = v. □ 
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Lemma 28. Suppose that I is the length of an optimal backward subpath of type 
I. In linear time, we can find an optimal backward subpath of type III or IV with 
length less than I or determine there is no such subpath. 

Proof. Similar to Lemma |2"5I the necessary and sufficient condition of type III 
shown in Lemma [27] can be determined by checking whether the max-flow in 
H£ is at least two or not. And the max-flow computations for all candidate 
pairs can be done in linear time. The optimal backward subpath of type IV can 
be computed similarly. □ 



4. Shortest detour path 

In this section we show an efficient algorithm for finding a shortest detour 
path. A shortest detour path contains exactly one outward subpath and has 
no backward subpath, in which an outward subpath is a path P such that 
E(P) C E — E(D), both endpoints of P are in V(D), and any of its internal 
vertex is not in V(D). Note that a simple si-path containing an edge not in D 
must have length strictly larger than d(s,t), or otherwise it should be entirely 
in D. Our goal is to efficiently find a minimum length st-path with an outward 
subpath. 

In this section T denotes an arbitrary shortest-path tree of G rooted at 
s and let F = T — E(D) denote the graph obtained by removing edges in 
E(T) fl E(D) from T. Apparently F is a forest consisting of subtrees of T and 
V(F) = V(T) = V. By the definition of D, any shortest path between two 
vertices in D must be included in E(D). For any v £ V(D), the path from s to 
uonT must be entirely within E{D) and therefore v must be a root of a subtree 
of F. Furthermore, the root of any subtree of F must be in V(D) because the 
edge between it and its parent is removed. 

Definition 7. For any vertex v £ V, let r v denote the root of the subtree of 
F which v belongs to. Let E denote the set of edges (x, y) such that (x, y) £ 
E - {E{T) U E{D)) and r x ^ r„. 
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Define 

| d s (x)+w(x,y)+d t (y) if (x,y) € E 
f(x,y) = < (2) 

I 00 otherwise. 

Note that, since G is undirected, both (x,y) and (y,a;) denote the same edge. 
But f(x,y) ^ f{y,x) in general. 

Lemma 29. ^Iny detour path P contains an edge in E. Furthermore, if {u, v) € 
£? is an edge on P, then f(u,v) ^ w(P). 

Proof. By definition P contains an outward subpath. Since the both end- 
points of this outward subpath are in V(D), they must be in different subtrees 
of F, and P must have an edge in E. The result f(u, v) ^ w(P) directly follows 
from definitions. □ 



For any vertex v, dt(v) ^ d(v, r v ) + dt(r v ) and the equality holds iff P[v, r v ] o 
<3[r„,£] is a shortest ui-path, in which P[v,r v ] is the ir^-path in T and Q[r Vl t] 
is an arbitrary shortest r„i-path in D. A vertex v is a dangler if dt(v) — 
d(v,r v ) + dt{r v ). By definition any vertex in V(D) is a dangler. 

Lemma 30. Ifv is not a dangler, there exists a detour path Q of length at most 

d s (v) +d t (v). 

Proof. Let P be the su-path in T and P' any shortest wt-path. Let q be the 
last vertex of P' intersecting P, possibly q = v. The path Q = P[s, q] o P'[q, t] 
is a simple si-path, and the length of Q is d s (q) + dt{q) ^ d s (v) + d t (v). We 
shall show d s (q) > d s (r v ). Then, by the definition of F, r q = r v and therefore 
q $l V(D). Consequently P is a simple path not entirely in D and thus a detour 
path. 

Suppose to the contrary that d s (q) Sj d s (r v ). Since dt(q) + d s (q) ^ d(s, t) = 
d s (r v ) + dt(r v ), we have dt{q) ^ dt(r v ) and furthermore dt{q) — d s {q) ^ d t {r v )~ 
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d s (r v ). Since q is on P', 

d t (v) = d(v,q) + d t (q) = (d s (v) -d s (q)) + d t (q) 
> d s (v) - d s (r v ) +d t (r v ) 
= d(v,r v ) + dt(r v ) 

which is a contradiction to that v is not a dangler. □ 



Lemma 31. Suppose that v is not a dangler and P is any shortest vt-path. For 
any u <E V(P), d s (u) > d s (r v ). 

Proof. Let u be a vertex with d s (u) ^ d s (r v ). We show that P cannot con- 
tain u. Since r v is on a shortest sv-path, d(v,u) ^ d(v,r v ), and therefore 
d(v,u) + d t (u) ^ d(v,r v ) + d t (u). Since d s (u) + dt(u) ^ d(s, t) = d s (r„) + iit(r l ,), 
by d s (u) ^J d s {r v ), we have dt(u) ^ rft(r. u ). Thus, d{v,u) + d t (u) ^ d{v,r v ) + 
dt(r v ). Since v is not a dangler, d(v,r v ) + dt(r v ) > dt(v), and therefore u is not 
on any shortest vi-path. D 



Lemma 32. If (x,y) minimizes function f and f(x,y) ^ oo, then there exists 
a simple st-path of length f(x,y) and with one edge in E. Such a path is a 
shortest detour path. 

Proof. Since an edge in E is not an edge in E(D), a simple path containing 
edge (x,y) € E must have length strictly larger than d(s,t). We only need to 
show the existence of such a simple path, and then it is a shortest detour path 
by Lemma l29l 

Let P x and P y be the shortest paths from s to x and y on T, respectively. 
Let Fbea shortest yi-path. If P x and Y are disjoint, P x o (x, y) o Y is a simple 
path and its length is clearly f(x,y). Otherwise let q be the last vertex of Y 
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Figure 5: (a) a; is a dangler while y is not. r x is a common vertex of a shortest sx-path and a 
shortest art-path; (b) The triangle inequalities when Y intersects P x (for Lemma l32ll ; (c) Case 
3 in the proof of Lemma 1321 when q is on P x [r x , x] . 



intersecting P x . By the triangle inequalities (see Fig. [5](b)), we have 

f(y,x) = d s (y) + w(y,x) + dt(x) 

< w(P x {s,q]oY[q,y]) + w(x,y)+w(P x {x,q]oY{q,t]) 

= w(P x ) +w(x,y) +w(Y) = d s (x) +w(x,y) + d t (y) = f(x,y). 

By the minimality of f(x,y), the equality must hold, i.e., 

d s (y) + d t (x) = d s (x) + d t (y). (3) 

We divide into three cases according to whether x and y are danglers. 

• Case 1: Both x and y are danglers. By Lemma [3J there are two disjoint 
paths Qi[s,r y ] and Q[[r Xl t\; or Q2[s,r x ] and Q^Vyjt] in D + . If Qi and 
Qi exist, the path Q\ o P y [r y , y] o (y, x) o P x [x, r x ] o Q' x is a simple path. 
Since X is a dangler, P x [x, r x ] o Q' x is a shortest xi-path and the length of 
the path is f(x,y). The other case that Qi and Q' 2 exist can be shown 
similarly. 
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• Case 2: Neither x nor y is a dangler. By Lemma OUl there exist two 
detour paths Pi and Pi such that w[P\) ^ d s (x) + d t (x) and w(P2) < 
d s (y) +d t (y). Then 

min{u;(P 1 ),u ) (P 2 )} < ( W (P) + w(P 2 ))/2 

< (l/2)(d.(*) + dt(a;)+d,(i/) + df(i/)) 

= d s (x)+d 4 (y)(byEq. ©) 

< f(x,y). 

• Case 3: Either x or y is a dangler. W.l.o.g. assume that £ is a dangler but 
y is not. First we show that d s {r y ) < d s (r x ) in this case. Recall that Y 
is a shortest j/i-path intersecting P x . If the intersection q is on P x [s,r x ], 
by Lemma 1311 d s (r y ) < d s (<7) ^ d s {r x ). Otherwise g is on P x [r x ,x] and 
therefore r x is on a shortest yt-path. As a result, d(j/, r x ) + dt(r x ) = dt(y). 
Since y is not a dangler, d(y,r x ) + d t (r x ) = dt(y) < d(y,r y ) + d t {r y ). 
Therefore 

d t (r x ) -dt(r v ) < d(y,r y ) - d(y,r x ). (4) 

Since r y is on a shortest sy-path, d s (r y ) +d(r y ,y) ^ d s (r x ) +d(r x , y), and 
equivalently 

d{r y ,y) - d{r Xl y) < d s (r x ) -d s (r y ). (5) 

Since both r^ and r a are in V(D), d s (r x ) + dt(r x ) — d(s,t) = d s (r y ) + 
dt{r y ), and we have 

d s (r a ) - d s (r x ) = d t (r x ) - d t {r y ) 

< d{y,r y )-d{y,r x ) (by Eq. ©) 
^ d,(r x ) - d a (r y ). (by Eq. ©) 

Therefore d s {r y ) < d s (r x ). Thus, any shortest r^t-path Q is disjoint to 
P y . The path P y o (y,x) o P x [^)?"x] ° Q is a simple si-path with length 
/(j/.a:) = f(x,y). 
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□ 



Theorem 33. For an undirected graph with nonnegative edge lengths, the short- 
est detour path problem can be solved in 0(m + n) time if d s (v) and dt(v) are 
given for all v. 

Proof. By Lemma [351 the length of a shortest detour path is the minimum 
value of function /. To compute {x,y) minimizing /, we first construct D and a 
shortest path tree T, and then find the edge set E. The minimum value of / can 
be found by checking both f(x, y) and f(y, x) for all edges (x, y) G E. The time 
complexity is linear if the distances d s (v) and dt(v) for all v are given. Once 
(x, y) is found, by the method in the proof of Lemma 1321 the corresponding 
path can be constructed in linear time. D 



5. Concluding remarks 

By Theorem [TJ we have the next corollary. 

Corollary 34. For undirected graphs with nonnegative edge lengths, if the sin- 
gle source shortest path problem can be solved in 0(t(m,n)) time, the next-to- 
shortest path problem can be solved in 0(t(m, n) + m + n) time. 

Important graph classes for which the single source shortest path problem can 
be solved in linear time include unweighted graphs (by BFS [3j), planar graphs 
[7[, and integral edge leng th graphs [171 ]. 

As pointed out in [12J LlSJj , it can be easily shown that the next-to-shortest 
problem is at least as hard as finding a shortest path between two vertices. 
When negative- weight edges are allowed, the next-to-shortest problem becomes 
NP-hard because it is polynomial-time reducible from the longest path problem 
by a similar reduction. An interesting problem is how to efficiently find the 
next-to-shortest paths for single source and multiple destinations. Another open 
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problem is the complexity of the version on directed graphs with positive- weight 
edges. 
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